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Abstract 

We revisit the standard treatment [Xu, Kharzeev, Satz and Wang, Phys. Rev. C 
53, 3051 (1996)] of J/ip suppression due to gluonic bombardment in an equilibrating 
quark-gluon plasma. Effects arising from gluon fugacity, relative g — ip flux, and ip 
meson formation time are correctly incorporated in the formulation of the gluon number 
density, velocity-weighted cross section, and the survival probability. Our new formulae 
are applied to numerically study the pattern of J/ip suppression in the central rapidity 
region at RHIC/LHC energies. The temperature and transverse momentum dependence 
of our graphs have noticeable differences from those of Xu et al. 

PACS numbers: 12.38M 
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1 INTRODUCTION 



Relativistic heavy ion collision experiments at CERN SPS/LHC and BNL/RHIC are 
believed to have led to a phase transition from the hadronic world into deconfined 
and/or chirally symmetric state of free quarks and gluons, the so called quark-gluon 
plasma (QGP) [l]-[6]. However, as of now, no conclusive evidence of QGP formation has 
been discerned. Among the most hotly debated, theoretically proposed signatures in this 
context are the erstwhile J/%1) suppression due to medium influence [6] and recent J/%1) 
enhancement via dynamical regeneration [7]. The well known mechanisms responsible 
for J/ ip suppression are summarized briefly in the Appendix for the sake of completeness. 

Attention in the sequel will be focused on the break-up of the J/ ip owing to bombard- 
ment with energetic gluons [8]. For this mechanism, Xu et al. [2] employed statistical 
mechanics coupled with phenomenological QCD to calculate the J/ip survival probability 
in a temporally evolving parton gas. The aim of the present paper is to extend/modify 
the work of Xu et al in the following multifold respects : 

(i) Gluon fugacity effect: For large momentum gluons which are responsible for 
the ip meson dissociation Xu et al. [2, Eq.(7)] replaced the gluon fugacity by unity 
in the denominator of the Bosc-Einstein distribution function. However, in the early 
stage of evolution the system may be quite far from chemical equilibrium implying that 
Xg need not be close to unity. In Sec. 2 below we derive a new formula for the gluon 
number density Ug valid for general Xg < 1. 

(ii) Relative flux effect : Xu et al. [2, Eq.(8)] were interested in the product F = 
Wreio" where the gluon--?/' break-up cross section a was written in the ip meson rest frame, 
but unfortunately their relative flux Vrei was evaluated in the fireball frame. In Sec. 3 
below we modify this procedure by treating the product F strictly in the ip rest frame. 

(iii) Format ion time effect : For computing their survival probability S{pt) Xu et 
al. [2, Eq.(14)] used an integration over having lower limit 0, where is the proper 
time measured in J/ip rest frame. This is very inconvenient because the gluon density 
ng{t) and the thermal- averaged cross section (vj-eicr) are natural functions of the usual 
time t in fireball rest frame. In Sec. 4 below we write a modified expression for S{pt) 
using t integration where formation times of the QGP as well as J/i/j are explicitly 
included. Of course, in our numerical results in Sections 2, 3, 4 the explicit velocity 
profiles of hydrodynamic flow are ignored. Finally, our main conclusions appear in Sec. 
5. 
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2 Number Density 



2.1 Preliminaries 

Assuming thermal equilibrium and working in the fireball rest frame let the symbol 
T denote the absolute temperature, K the gluon four momentum, 16 the spin-colour 
degeneracy factor, A^, < 1 the gluon fugacity, and / = f{K^,T,\g) the one-body gluon 
distribution function. Then the gluon number density Ug is obtained from 



2.2 Xu Procedure 



For near chemical equilibration Xu et al. [2, Eq.(7)] employed an approximate, factorized 
Bose-Einstein distribution 



- = A. E e-"""/" (2) 



^ n=l 

which led them to a number density depending on the fugacity linearly through 



< = f,T%a^) (3) 



2.3 Our Proposal 

In order to tackle the possibility of gluon chemical non-equilibration we use the full 
Bose-Einstein form 

\ oo 

fOur _ _ \n -nK°/T (a\ 

■I ~ f,KO/T _ \ 2^ ^9^ 

^9 n=l 

which guides us to a number density containing the fugacity in a power series via 



n 



Our 



16 , ^ A" 



= zi^' E -I (5) 



Remembering that l/rr" type scries converges rapidly with n, numerical comparison of 
Eqs. (3, 5) is easily done via the ratio 



which, however, requires the knowledge of at various times. Of course, the algebraic 
reason for the inequality rf^^ < n^^ is the fact that the distribution function /^^"^ < 
as long as Xg < 1. In other words, the effect of correct fugacity (before chemical 
equilibration) is to reduce the number density of gluons below the value of Xu et al. 



2.4 Initial Conditions 

It is well recognized that the scenario resulting from relativistic heavy-ion coUisions is 
rapidly time-dependent. Quick scattering among the partons drives the QGP to thermal 
equilibrium in the fireball rest frame within a time = Tq ~ 1/A ~ 0.5fm/c where the 
suffix i stands for "initial" and A is the QCD energy scale. The initial conditions 
predicted by HIJING Monte Carlo simulation are summarized in Table 1. There gluon 
densities computed via Xu procedure (Eq.3) and our proposal (Eq. 5) are also listed. 
Clearly the relative difference between n^^ and n°"'' is of the order of 1/8 ~ 12% which 
is significant. 

Table 1: Initial values for the temperature, time, fugacities etc. at RHIC(l), LHC(l) 
only [5] 





T (GeV) 


ti = To (fm) 


A. 


\ 


n^"(fm)"^ 


nO^"-(fm)-^ 


RHIC(l) 


0.55 


0.70 


0.05 


0.008 


2.11 


1.76 


LHC(l) 


0.82 


0.5 


0.124 


0.02 


17.34 


14.66 



2.5 Temporal Evolution 

The thermally equilibrated QGP produced at the instant ti — tq undergoes rapid ex- 
pansion (accompanied with cooling) while partonic reactions tend to drive the plasma 
towards chemical equilibrium. In Bjorken's boost-invariant longitudinal expansion sce- 
nario the fugacities and temperature are known [16] to evolve through the following 
master rate equations : 

^ + 3| + ^ = i?3(l-A,)-2i?2 (l- 

Xq f 1 
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b \ 

\ + —\] '^^^ = const (7) 

Here r is the medium proper time, Ag the quark fugacity, Nf the number of flavours, 
and remaining symbols are defined by 

ai = 16C(3)/7r^ as = Sir^/lb 

bi = 9C{3)Nf/7r\ 62 = 77r'^Nf/20 (8) 

Their solutions on the computer yield the functions T{t), Xg{t), ng{t) in terms of the 
fireball time t. The lifetime (or freeze-out time) tufe of the plasma is the instant when 
the temperature drops to T{t^.fJ — 200 MeV, say 



3 Flux-weighted Rate 
3.1 Preliminaries 

Next, the question of applying statistical mechanics to gluonic break-up of the J/ip 
becomes relevant. In the fireball frame consider a ip meson of mass m^, four momentum 
p^, three velocity — p-tp/p^, and dilation factor 7^ = p^/m^. If q is the gluon four 
momentum measured in ^jJ meson rest frame then by Lorentz transformations 

= 7i/'(9° + q) = 7v5°(l+ I ^V- I cos9q^) 
£K = {Kyq°) d?q (9) 

where 9q^ is the angle between q and unit vectors. 

The invariant quantum mechanical dissociation rate for g — ip collision can be written 
compactly as 

r = v^^, a (10) 

where Vrei is the relative flux and a the cross section written in any chosen frame. Its 
thermal average over gluon momentum in fireball frame reads 



^ ' UgJ (27r)' q^ ^ ^ ' 

with / being the distribution function already encountered in Eq.(l). 
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3.2 Xu Procedure 



Xu et al. [2, Eq.(8)] worked with their relative flux v^f = qViK'^-ii,) in the fireball 
frame, but unfortunately the cross section a^^^ was in meson rest frame based on the 
standard QCD value [20] 

"n^ = BiQ" - iflQ"" (12) 

«°4'^=f(fyV-W. (13) 

where is the J jil) binding energy and mc the charm quark mass. Insertion into Eq. 
(11) led them to 

where the approximate given by Eq.(2) has been recalled. The simple angular 
integration over d cos Qq^ can be done by taking the polar axis along to yield 



T'i-rs/i \ ^^ . I T7XU 

"5 n 



00 -1 

E ^ / dQ^ gV,,,, (e-^"-^° - e-^"^«°) (16) 



^ li, I 

Here the following abbreviations have been introduced : 

=1 i;^ I a , = Cn ± = (1± | |) (17) 

For J/i/j produced in the central rapidity region, Xu et al. have drawn elaborate curves 
showing the dependence of (P-^") on T and pt- 



3.3 Our Proposal 

We set up our F entirely in the ip meson rest frame where v^^f^ — c — 1. Thereby Eq. 
(11) becomes 

(r°"0 = / 7f^7^ (1+ I I cos e,^) a,,, f: A^e-^"/^ (18) 



riOur 

"-g 
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where the exact /°"'' given by Eq.(4) has been recalled. The slightly complicated angular 
integration over d cos 9q^ can be performed by choosing the polar axis along yielding 



7r^\v^\ nO""- — 1 n Ji 



1- \v^\ { 



1 - -)] e'^^^° 



l+\v^\{l + —) \ e 

K Pn J 



(19) 



where 



Pn = DnQ^ = D„q^/e^ , Io{pn) = (sinhp„)/p„ 
h{Pn) = (coshp„)/p„ - (sinhp„)/p^ 



(20) 



Clearly the dependence of on Xg and 7^ is more involved than that of (F ) given 

by Eqs. (15, 16). 



3.4 Numerical Work 

The initial thermally-averaged rates (F-^") (Eq.l6) and (^F*-*"') (Eq.19) are depicted in 
Figs. 1, 3 and Figs. 2, 4, respectively. The physics of dependence of the peak on the 
temperature and transverse momentum has already been discussed in Ref. [2]. Here we 
wish to focus attention only on the striking similarity between Figs. 1 and 2 inspite of 
the different fugacities and fluxes employed. For this purpose, we first go back to the 
Lorentz transformation (9) and observe that the Xu et al relative flux receives dominant 
contribution from the antiparallel (cos^g^ = —1) configuration. Indeed, then 



7^X0 7^(1- I I) 



1+ I 



(21) 



Now let us consider the ratio 



(FO" 



n 



Xu 



phase space integral of a^^^^^ 



phase space integral of c cTj, 



(22) 



Due to fugacity effect the number density ratio in Eq.(22) is somewhat smaller than 
unity as already mentioned in Sec.2. On the other hand, since in near antiparallel 
configuration, the relative fiux > c, hence the ratio of the phase space integrals is 
somewhat larger than unity. These two effects tend to partially compensate each other 
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Figure 1: The thermal-averaged gluon-J/ip dissociation cross section {v^^^a) as a function 
of temperature at different transverse momenta as done by Xu et al. [2, Eq.(16)]. The 
initial gluon fugacity is given in Table 1 at RHIC energy. 



in Eq.(22) so that the relative difference between the curves of Figs. 1 and 2 is not more 
than about 5 — 6%. However, the influence of 1+ | \ becomes more pronounced at 
high transverse momentum, causing noticeable difference between the curves of Figs. 3, 
4. 



4 Survival Probability 
4.1 Preliminaries 

Consider a cylindrical coordinate system in the fireball frame where the ip meson was 
created at the time-space point {ti,r^, (plp) with transverse velocity if^. The plasma is 
supposed to be contained within a cylinder of radius R, expanding longitudinally till the 
end of its lifetime t^.^^. The ip meson's trajectory will hit the said cylinder after covering 
a distance o?^^ in the time interval such that 

dm = -r^ cos 0^ + yi?2 _ 2 ^.^2 ^7 

tm = d^i/ \ ^\ > (23) 
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Velocity-weighted cross section (Our Model) 
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Figure 2: Same as in Fig.l, but these curves are obtained by our Eq.(19). 

the full temporal range of interest is obviously 

ti <t<tii ; tii = min {U + tm , tM^) 



(24) 



The corresponding survival probability of J/ip averaged over its initial production con- 
figuration extending over the transverse area A becomes 



S{pt) = d% {r' - rf ) e-^/ d% {r' - rf ) 
W = (" dtugit) {T(t)) 

Jtr 



(25) 



4.2 Xu Procedure 

Xu et al. [2, Eq.(14)] did not take into account the formation time of the coulombic 
bound state, i.e., they chose the instant of production as 



tf^ — ti — To 



(26) 



Also, they seem to have used as integration variable the proper time = {t — ti) /^^ 
measured in meson rest frame. This procedure is inconvenient since the gluon number 
density n^" was best known in the fireball frame. 
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Velocity-weighted cross section (Xu et al) 
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Figure 3: The thermal-averaged gluon- J/ip dissociation cross section {v^^^a) as a function 
of transverse momentum at different temperatures as done by Xu et al. [2] . 

4.3 Our Proposal 

We do take into account the formation time ~ 0.89fm/c of the bound state in the cc 
barycentric frame. Remembering the dilation factor 7^, we choose 



and retain the fireball time t for integration in Eq.(25). 
4.4 Numerical Work 

In Figs. 5, 6 the J/ip survival probabihty has been plotted as a function of the transverse 
momentum based on the general formula (25). The solid curve denotes our result using 
Eq.(27) while the dashed curve is that of Xu et al employing Eq.(26). Clearly, the 
J/ip^s survival chance is much more (i.e.. their suppression is substantially less) in our 
case compared to Xu et al's. Its reason can be understood by examining the integral 
appearing in Eq.(25) viz. 



Our 



ti + lijTp 



(27) 



W = 



/ dt [phase space integral of Vj-eiCJ over / at time t\ 
ti 



(28) 
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Velocity-weighted cross section (Our Model) 
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Figure 4: same as in Fig. 3 except the results are obtained by us. 



First, we recall from Sec. 2 that /*^"'' < /"^^ due to the fugacity effect. Next, we know 
from Eq.(21) that v^ei = c < v^f due to the flux effect. Finally, Eqs.(26, 27) tell that the 
time interval available for dissociation tu — tf^^ < tu — tf^ due to the formation time 
effect. These three mechanisms operate cooperatively to make VF*-*""^ < 1^-^" resulting in 
less suppression. 



5 Conclusions 

(i) In this paper we have extended the work of Xu et al [2] concerning the gluonic break- 
up of the J/ip^s created in an equilibrating QGP. Our theoretical formulae on number 
density (Eq. 5), flux- weighted cross section (Eq.l9), and survival probability (Eq. 27) 
are new. 

(ii) Our numerical results are also significant as compared to those of Xu et al. Since 
gluon fugacity is less than unity before chemical equilibration, hence our number density 
n^"'^(i) of hard gluons (which are primarily instrumental in dissociating the J/ip^s) is 
lower as shown in Table 1. 

(iii) Next, since our g — ip relative flux in meson rest frame is only v^J^'^ = 1 (and 
not H- I -y^ I of the fireball frame) hence our thermally- averaged rate (r'-*"''(t)) is also 
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Figure 5: The survival probability of J/ if) in an equilibrating parton plasma at RHIC 
energy with initial conditions given in Table 1 [2]. The solid curve is our result, while 
the dashed curve is the result obtained by Xu et al. [2] . 

smaller as depicted in Fig.4. 

(iv) Since we properly take into account the production time of the J/ip^s, hence the 
temporal span available for their break-up becomes shorter. These three effects act in 
a cooperative manner to reduce substantially the amount of J/ip suppression (i.e. to 
increase noticeably their survival chance S'^^'^{pt) as demonstrated by Figs. 5, 6. 

(v) Apart from possible cc recombination [7] another important effect not considered 
in the present paper is the transverse hydrodynamic expansion of the QGP. Mathemat- 
ically such expansion demands that the gluon statistical mechanics must be done in 
a local comoving frame, while physically the temperature will drop more quickly with 
time. This highly nontrivial problem is under investigation at present and its results 
will be published in a future communication. 
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Figure 6: Same as in Fig. 5 but it is at LHC energy. 
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Appendix 



J/ifj Suppression Mechanisms Summarized 

In relativistic heavy ion collision the heavy quark-antiquark pairs (leading potentially 
to J I'll) mesons) are produced on a very short time scale ~ l/2mc — 10^^"^ sec with rric 
being mass of the charmed quark. The pair develops into the physical resonance over 
a formation time ~ 0.89fm/c in its own rest frame. This J/V' traverses the deconfined 
plasma together with the hadronic matter before leaving the interaction region to decay 
into a dimuon which is finally detected. However, this chain of events can be prevented 
via any of the following mechanisms. 

Even before the cc bound state is created it may be absorbed by the nucleons stream- 
ing past (Glauber/normal absorption [6]). Or, by the time the resonance is formed the 
Debye screening of the colour forces in the plasma may be sufficient to kill it [9] . Or, an 
energetic parton could hit and dissociate the J/ip [8]. Or, the Brownian motion of the 
J/ip through the medium could cause its sufficient swclling/ionization [3]. The extent 
of suppression will be decided by a competition between the J/ip momentum and the 
rate of hydrodynamic expansion (with associated cooling) of the plasma [11]. Of course, 
the entire above picture will be substantially modified if the J/^'s are regenerated via 
cc recombination [7] . 
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